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THE CRITICAL EXPONENTS OF THE QCD (TRI)CRITICAL ENDPOINT 
WITHIN EXACTLY SOLVABLE MODELS 
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Bogolyubov Institute for Theoretical Physics, Kiev, Ukraine 

The critical indices a' , /?, 7' and S of the Quark Gluon Bags with Surface Tension Model with the tricritical and 
critical endpoint are calculated as functions of the usual parameters of this model and two newly introduced 
^^ parameters (indices) . The critical indices are compared with that ones of other models. The universality class 

I of the present model with respect to values of the model parameters is discussed. The scaling relations for the 

f^ found critical exponents are verified and it is demonstrated that for the standard definition of the index a' some 

CN of them are not fulfilled in general case. Although it is shown that the specially defined index a'^ recovers the 

1—t scaling relations, another possibility, an existence of the non-Fisher universality classes, is also discussed. 
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1 Introduction 



,_i Investigation of the properties of strongly interacting matter equation of state has become a focal point of 

1-^ modern nuclear physics of high energies. The low energy scan programs performed nowadays at CERN SPS 

Mh and BNL RHIC are aimed at the discovery of the (tri)critical endpoint of the quantum chromodynamics (QCD) 

O . phase diagram. Despite many theoretical efforts neither an exact location nor the properties of the QCD 

D (tri)critical endpoint are well known [T]. Therefore, the thorough theoretical investigation of the QCD endpoint 

;^ , properties are required in order to clarify whether this endpoint is critical or tricritical. 

This work is devoted to calculation of the critical exponents of the QCD endpoint of both critical (CEP) and 
^-H tricritical (triCEP) types. Unfortunately, such a task cannot be solved within the QCD itself. Moreover, even 

<^ the possibilities of the lattice QCD are nowadays very limited in this respect. Then, unavoidably, one has to use 

^-j some models of the (tri)CEP. The most popular models of this kind are the quark-meson model [HE] and the 

^k extended Nambu-Jona-Lazinio model 4J . Despite a great popularity both of these models contain a significant 

^-H weakness, they are mean-field ones and hence there is no reason to expect that their critical exponents could 

I — I differ from that ones of the Van der Waals model equation of state [51. Therefore, to study the properties of 

f— ^ non-classical endpoints one has to investigate the non-mean-field models. 

T-H The gas of bags model [6 is an example of such a non-mean-field model. Despite a significant success of 

^"^ this model in describing of the deconfinement phase transition, it can not generate the (tri)CEP in a natural 

^ way. The problem was solved in the Quark Gluon Bags with Surface Tension Model (QGBSTM) [7l[8]. This 

exactly solvable model accounts for the surface effects which play a decisive role in the critical phenomena and 
employes the same mechanism of the (tri)CEP generation which is typical for the liquid-gas phase transition 
(PT) and which is also used in the Fisher droplet model (FDM) [5] and in the statistical multifragmentation 
model (SMM) [TUl E]: the endpoint of the 1-st order PT appears due to vanishing of the surface tension 
coefficient at this point which leads to the indistinguishability between the liquid and gas phases. However, the 
surface tension coefficient in the QGBSTM has the region of negative values, which is a principally different 
feature of this model compared to the FDM, SMM and all other statistical models of the liquid-gas PT. Note 
that just this feature provides an existence of the cross-over at small values of the baryonic chemical potential 
fi in the QGBSTM with triCEP [7j and with CEP [TT], and also it generates an additional PT in the model 
with triCEP at large values of fj,. Therefore, it is very important and interesting to study the critical indices of 
such a novel statistical model as the QGBSTM, to determine its class of universality and to examine how the 
latter is related to that ones of the FDM and SMM. 

The paper is organized as follows. A description of the QGBSTM is given in Section [2] In Section [3] the 
critical exponents of the CEP and triCEP are calculated. This section is also devoted to the analysis of scaling 
relations between the found critical exponents. Conclusions are given in Section |4] 

2 Quark Gluon Bags with Surface Tension Model 

An exact solution of the QGBSTM was found in [7]. The relevant degrees of freedom in this model are the quark 
gluon plasma (QGP) bags and hadrons. The attraction between them is accounted like in the original statistical 
bootstrap model [12] via many sorts of the constituents, while the repulsion between them is introduced a la 
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Van der Waals equation of state Oil]. An essential element of the QGBSTM is the T and /i dependence of its 
surface tension coefficient T'E{T,ii) (here E(r, /z) is the reduced surface tension coefficient). Let us denote the 
nil line of the reduced surface tension coefficient in the T — ^ plane as T^{fj,), i.e. 2(7^, /i) — 0. Note that for a 
given fi the surface tension is negative (positive) for T above (below) Ts(/u) line. Here it is appropriate to say a 
few words about the negative values of S which is a distinctive feature of QGBSTM compared to other models. 
There is nothing wrong or unphysical with the negative values of surface tension coefficient, since TEw" is the 
surface free energy of the bag of mean volume v and, hence, as any free energy, it contains the energy part Csurf 
and the entropy part Ssurf multiplied by temperature T [5]. Therefore, at low temperatures the energy part 
dominates and surface free energy is positive, whereas at high temperatures the number of bag configurations 
with large surface drastically increases and it exceeds the Boltzmann suppression and, hence, the surface free 
energy becomes negative since Ssurf > '^^ ■ Such a behavior of the surface free energy can be derived within 
the exactly solvable model of surface deformations known as Hills and Dales Model [13]. 

In the grand canonical ensemble the pressure of QGP and hadronic phase are, respectively, given by 

PQ{T,fi) = Tsq{T,^), (1) 

PH{T,fi) = TsHiT,fi)=T[FH{sH,T,fi)+uiT,fi)Ir{As,E)], (2) 

n °° 

FH{sH,T,fi) = ^5,e^-'^^-^-0(T,m,) and /.(As, S) = I [|^e-^™-^'^^ (3) 
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where As = sh(T, ii) — sq(T, fi). The particle density of a hadron of mass rrij, baryonic charge bj, eigenvolume 
and degeneracy gj is denoted as (t>j{T,mj) = ^^ Jp^dpe t . Here, as in [7], it is assumed that 
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the functions u(T,ii) and sq(T,ii), which are the parameters of the present model, and their first and second 
derivatives with respect to T and /i are finite everywhere at the T — fi plane. In the continuous part of the 
spectrum of bags denoted as u/t(As, E) the surface of a QGP bag of volume v is parameterized by the term 
v'^. Usually one chooses k = | in 3-dimensional case (or n = ^^ for the dimension d), but in what follows it 
is regarded as free parameter of the range < k < 1 . 

The system pressure, that corresponds to a dominant phase, is given by the largest value between pQ and 
Ph for each set of T and /i. As usual, the deconfinement FT occurs when pressure of QGP gets equal to that 
one of the hadron gas, which is nothing else as the Gibbs criterion. The sketch of T — /x phase diagram with 
triCEP is shown in the left panel of Fig. [l] Suppose that the necessary conditions for the deconfinement FT 
outlined in [7] are satisfied and its transition temperature is given by the function Tc{fi) for /i > ficep- The 
necessary condition of the triCEP occurrence is that at the nil line of the surface tension coefficient there exists 
the surface tension induced FT of 2-nd (or higher order) f7| for /i > fjLcep and Tc{ii) < TY,{ii) for these fj, values. 
The both of these inequalities become equalities only at triCEF. Moreover, at triCEF the phase coexistence 
curve Tc{iJ,) is a tangent (not intersecting!) line to the nil line of the surface tension coefficient T^{fi). For 
fi < iicep the deconfinement FT degenerates into a cross-over since in this region E < |7j and, hence, the 
system pressure is defined by a solution of Eq. (l2|. It is necessary to stress here that QGF exists as real quark 
gluon liquid (QGLiquid) with the pressure (II]) only for Tc{ii) <T< T^ifJ.) and /i > /icep, outside this region the 
system is defined by Eq. (pi). Now let us discuss phase diagram with CEF (see the right panel in Fig. fTl). The 
only difference with the case of triCEF is that Tc{fi) and Tj^ip) lines coincide exactly for fj, > ficep- Therefore, 
QGLiquid exists only at the FT line. 

An actual parameterization of the reduced surface tension coefficient E(T, /i) is taken from [71ll4j: 



E(T,/x) 



T 



rs(M)-T^ 



T^{^^) 



Cn /„ „ dT, 

TV dfi 



E (T, /i) = ^— ■ ( Tcep - r + ——{flcep - fJ-) 



sign (7s ip) - T) for triCEF , (4) 



T^{^A 



for CEF . (5) 



In what follows the coefficient ctq is assumed to be a positive constant, i.e. (Tq = const > 0, but it is easy to 
show that the obtained results hold, if ctq > is a smooth function of T and /x. 

In the vicinity of (tri)CEF the behavior of both the deconfinement FT curve and the nil surface tension 
coefficient line in the n — T plane is parameterized via a single parameter ^^ > 0: 

Tcep - Tt,{I-1.) ^ (a* - ^J^cep)^ and Tcep - Tc{fJi) -- (/i ~ fJ-cep)^ , (6) 

since, as discussed above, Tc{fJ.) and T^{p) lines are tangent to each other at the endpoint in the model with 
triCEF and coincide in the model with CEF. This is a new index which was not considered both in the FDM 
and SMM and as we show below, it is responsible for a new universality class compared to other exactly solvable 
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Figure 1. A schematic phase diagrams with triCEP (left panel) and CEP (right panel) in the fJ. — T plane. 
The dashed curve indicates the nil line of the surface tension coefficient Ts(/u), below (above) which S is positive 
(negative). The deconfinement FT line Tc{fJ.) is shown by the full curve for /i > ficep- At the phase diagram 
with CEP Tc{fJ,) and Ts(m) lines coincide exactly for /x > Hcep- A cross-over (the short dashed curve) takes 
place along the line T^i^i) for ^ < fXcep- The cross-over and PT regions are separated by (tri)CEP (filled circle). 



Since the entropy density and the baryonic density are, respectively, defined as T and fi partial derivatives 

for 



of the corresponding pressure, then using (1| and (2) the Clapeyron-Clausius equation -J^ = " _ '^ 

the model with triCEP can be explicitly rewritten as 



T = Tr 



dT 



At 

'a7- 



.(0,2) 






(7) 



T=Tc 



where the following notation is used A, 



- dFn 



ai 



dsQ 
di 



. ds 



1) for i e {T,fi}. Note that for the model 



with CEP the term ^u/^_k(0, S) does not appear in (|7l). Let us parameterize the behavior of the numerator 
and denominator in (|7| at the (tri)CEP vicinity as 



An 



\T=Ta 



{Tcep - T,{fl)f+ (' 



and A 



t^ I T=Ta 



(Tcep - Tc{fJ,)y 



(8) 



where x ^ max(0, 1 — tr ) denotes another new index. The latter inequality follows from the fact that the 
integral /,- and functions F^, u, sq together with their derivatives are finite for any finite values of T and fi. 
An introduction of the index x is quite general and could be done for any model with the PT of the liquid-gas 
type, since it unavoidably appears from the Clapeyron-Clausius equation, which is a direct consequence of the 
Gibbs criterion for phase equilibrium. To our best knowledge, this index was never used for the calculation of 
critical exponents. 

3 Critical exponents of the QGBSTM 

The standard set of critical exponents a' , /3 and 7 [SJ[TS] describes the T-depcndence of the system near (tri)CEP: 



Ap 
AKt 



|i|-" , for t<0 and p 
\tf, for t < 0, 
|i|-'^', for t<0. 



Peep: 



(9) 
(10) 

(11) 



where Ap — 



— ph)t=t^ defines the order parameter, Cp = — (f^)p denotes the specific heat at the critical 



density and AKt = {Kt — K!f)T=Ta is the discontinuity in the isothermal compressibility Kt = l(g^)T across 



the PT line, the variable t is the reduced temperature t = 
index 5 [HI [5] (hereafter the tilde indicates that T ~ Tcep 



T-Ta, 



The critical isotherm shape is given by the 



Peep - p ~ (pcep - pf for t = 0. 



(12) 



In this work we want to present the results of our calculations which are exhaustively explained in |16| . 
Not going into details here we give the standard set of critical exponents (see Table [I]) and discuss the related 
physics issues. Since in some aspects the QGBSTM is similar to the FDM and SMM it is interesting to compare 
its critical exponents with that ones of the FDM 9J and SMM I?! . The QGBSTM with CEP generates quite 
unique set of the critical exponents which are independent of r. Therefore, this model belongs to different 
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Table 1. Critical exponents of QGBSTM. For the model with CEP the notation /3±= C^+C'^-gr is introduced. 
The index a'^ describes the specific heat difference AC = [Cp^ — CpQ)T=Ta ~ \i['°'° foi' two phases [T3] . 

two regimes switched by parameter x- Since the FDM and SMM implicitly treat the parameter x = 0, then it is 
most natural to compare their critical exponents with the QGBSTM results just for this case. In the QGBSTM 
with triCEP there is a regime, when its index 5\x=o = §Ef matches the SMM result [T7|. Moreover, it is easy 
to see that in this regime all other indices the QGBSTM and SMM coincide for f "^ < 1. The spectrum of values 
of the QGBSTM critical indices is more rich than the corresponding spectra of the FDM and SMM since this 
model contains two new indices ^^ and x- 

The explicit expressions for the QGBSTM critical exponents allow us to examine scaling relations between 
them. For the model with triCEP Fisher inequality (a' + 2/3 + 7' > 2) and Griffiths inequality {a' + l3{l + 6) > 2) 
are not fulfilled in a general case, whereas the Liberman inequality {j' + /3{1 — S) > 0) is always obeyed \TE\ . 
For the model with CEP the situation is similar. In order to 'save' the scaling inequalities it was suggested to 
replace the index a' by a'^ [TS] . This hypothesis resolves the problem of scaling inequalities only for the triCEP 
model. Therefore, it is quite possible that the non-Fisher universality classes exist. 

4 Conclusions 

The critical indices of the QGBSTM are found in terms of the model parameters. Two parameters ^-'" > and 
X > max(0, 1 — -Ft) are newly introduced and consequently the spectrum of the values of critical exponents of 
the present model is more rich compared to other models. It is shown there is a regime when QGBSTM with 
triCEP reproduces the critical exponents of the SMM with triCEP, whereas the critical exponents of the FDM 
are never reproduced by QGBSTM. Thus, for x = QGBSTM and SMM fall into the same universality class. 
The direct calculations show that for the standard definition of the critical index a' (found along the critical 
isochore) the Fisher and Griffiths scaling inequalities are not always fulfilled, whereas the Liberman inequality 
is obeyed for any values of the model parameters. According to Fisher hypothesis the scaling relations for the 
index a^ are verified. The possibility of the non-Fisher universality classes existence is discussed. 
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